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The fundamental quantities of potential scattering theory are generalized to accommodate long- 
range interactions. New definitions for the scattering amplitude and wave operators valid for arbi- 
trary interactions including potentials with a Coulomb tail are presented. It is shown that for the 
Coulomb potential the generalized amplitude gives the physical on-shell amplitude without recourse 
to a renormalization procedure. 
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I. INTRODUCTION 

The two-body scattering problem is a central subject of 
quantum mechanics. It is well known that conventional 
quantum collision theory is formally valid only when the 
particles interact via short-range potentials (see, e.g. Q). 
In the time-dependent formulation, formal scattering the- 
ory can be made to include Coulomb long-range poten- 
tials by choosing appropriately modified time evolution 
operators 0, Q • This is equivalent to choosing various 
forms of renormalization methods 0, El El SIS 13 m the 
time-independent formulation. 

Though the renormalization theories lead to the cor- 
rect cross sections for the two-body problem, the re- 
sults from these procedures cannot be regarded as com- 
pletely satisfactory. For instance, in screening-based 
renormalization methods 0, different ways of shield- 
ing lead to different asymptotic forms for the scattering 
wave function. Generally, these asymptotic forms dif- 
fer from the physical one obtained from the solution of 
the Schrodinger equation [ljj. The weakest point about 
these methods, however, is that they give rise to a scat- 
tering amplitude that does not exist on the energy shell. 
This is because the amplitude obtained in these methods 
has complex factors which are divergent on the energy 
shell 0[ll|,mEE10- These factors, usually contain- 
ing branch point singularities, must be removed (renor- 
malized) before approaching the on-shell point. Further- 
more, the renormalization factors depend on the way the 
limits are taken when the on-shell point is approached. 
Thus, the ad-hoc renormalization procedure is based on 
prior knowledge of the exact answer to compare with and 
has no ab initio theoretical justification. These issues are 
discussed in detail in the comprehensive coverage of the 
subject given by van Haeringen |10| . 

The motivation for the present work is to demon- 



strate that there is a practical approach to the two- 
body collision problem with a Coulomb-like potential 
that does not lead to the formal difficulties described 
above. Our approach is based on a representation of the 
scattering amplitude written in a divergence- free surface- 
integral form which is ideally suited for practical calcu- 
lations. We build on a recent formalism which has im- 
proved our understanding of the three-body scattering 
processes 0, . In this work we consider local poten- 
tials. It is not difficult to show however that the present 
results hold for nonlocal potentials as well. 

In Section II we present a formal solution of the 
Schrodinger equation which satisfies all necessary condi- 
tions imposed by the long-range nature of the Coulomb 
interaction. We introduce new well-defined forms of the 
scattering amplitude and wave operators for two-body 
systems valid for arbitrary interactions. The new defini- 
tions of the scattering amplitude not only cover arbitrary 
potentials but also directly give the physical result. The 
relationship with conventional formulations is discussed 
in Section III. Section IV contains some concluding re- 
marks including a brief discussion of the utility of the 
two-body formalism for three-body Coulomb scattering 
problems above the breakup threshold. 



II. FORMALISM FOR LONG-RANGE 
POTENTIALS 

We consider scattering in a system of two particles 1 
and 2 interacting via an arbitrary spherically-symmetric 
potential V with the Coulomb long-range tail. Through- 
out the paper we use such units that % = 1. A scattering 
state of this system is the solution to the Schrodinger 
equation 



(E-H)1>£(r) = 0, 



(1) 
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where H = Hq + V is the total two-body Hamiltonian 
of the system, Hq = —A r /2[i is the free Hamiltonian, 
E = fc 2 /2/i is the energy of the system, r is the relative 
coordinate of the particles 1 and 2 and k is their relative 
momentum, /j, is their reduced mass. To be more specific 



2 



we can assume that interaction V consists of some short- 
range part V s and the Coulomb potential V c — ziz 2 /r, 
where z\ and zi are the charges of the particles. 

From all possible solutions to Eq. we should choose 
the one satisfying the asymptotic boundary condition 
corresponding to the physical scattering picture. When 
the potential has the Coulomb tail the scattering wave 
function tf)t(r) asymptotically behaves, in the lead- 
ing order, like the Coulomb-modified plane wave and a 
Coulomb-modified outgoing spherical wave 



1#(r) 



e ik-r+i~i ln(fer-fe-r) ^ _|_ 0(1 M] 
gikr — i'y ln(2fcr) 



-f(k ■ r) 



[l + 0(l/r)],(2) 



where 7 = z\Z2[i/k and / is the scattering amplitude. 
The second suitable solution ip^ (r) asymptotically be- 
haves, in the leading order, like the Coulomb-modified 
plane wave and a Coulomb-modified incoming spherical 
wave 

^ki r ) r ^°° e lk - r -^ x ^ kr+k - r \l + 0{l/r)] 



+f*(-k-r)- 



-ikr-\-i-f ln(2fcr) 



[l + 0(l/r)]« 



(3) 



Note that k ■ r ^ ±1, respectively for J2J an d ©■ 
However, as we will see below, in the asymptotic sense 
for which Eqs. (0) and J2J are written these for- 
ward/backward singularities have a 8- function nature 
and are, therefore, integrable. We will also see that seem- 
ingly different forms of the Coulomb- modified plane wave 
in Eqs. J2J and J2J asymptotically are essentially the 
same function. 

We can separate ip^ into the so-called 'incident' and 
the 'scattered' parts according to 



(4) 



where cjx^ and Xk asymptotically behave like the first and 
the second terms of Eqs. and ©, respectively. The 
'unscattered' wave incident at infinitely large distances is 
given by H3 



<t>i{r) 



[(W)«] 2 (Tikr + ik ■ r)- n /n\, (5) 



where (z) n = z(z + 1) • • • [z + n — 1). As to the un- 
known scattered wave xi-> the f° rm of Eqs. J2J) and © 
suggest that its leading order term in the asymptotic re- 
gion already contains all the scattering information we 
want. The next order terms simply repeat this informa- 
tion. Therefore, all we need for extracting the scattering 
amplitude is the leading-order asymptotic term of the 
scattered wave xt- 

Let us denote <jyi , (t") the first term of the incident 



wave 



k • 



<t> k ( r ) 



(6) 



and by 4% \ r ) the second term etc.: 



0L° )± (r) + ^ 1)± (r) + ^ a)± (r) + ..-. (7) 



Then Eq. can be written in the form 

(E-H)fy£(r)-ct> k 0)± (r)\ = (H - E)^ ]± (r).i s ) 
If we introduce Green's function according to 

(E-H)G(r,r',E) = 5(r-r'), 
and apply it onto both sides of Eq. (JSJ we have 

dr'G(r,r';E±ie)(E-H) (r') - 4 k 0)± (r') 

dr'G(r, r'\ E ± ie)(H - £)<^ 0)± (r'),(10) 



(9) 



where e is a small positive parameter and limit as e — > 
is assumed ^8|. We used an arrow on the differential 
Hamiltonian operator to show the direction in which it 
acts. We emphasize here a subtle point that the oper- 
ator G(r,r',E)(E - H) in Eq. ijlfjj) does not act like a 
8- function. In other words, though 

G(r, r', E)(E -H) = (E- H)G(r, r', E) = 8(r - r% 

(11) 

however, in general 

G(r, r', E)(E -H)^ 8(r - r'). 



(12) 



The reason is that the operator G(r, r', E)(E — H) pro- 
duces an integral that has a surface- integral component. 
In order for the action of the operator G(r, r', E)(E — H) 
to be equal to that of G(r,r',E)(E — H) this surface- 
integral component should be zero. A similar problem 
has been discussed by Glockle [19| in relation to the 
three-body Lippmann-Schwinger equations. Using an 
operator similar to ours but with the Hamiltonian H 
and the Green's function G replaced by some channel 
Hamiltonian and corresponding Green's function Glockle 
showed that surface integrals of the same nature only 
disappear if the operator acts on a function which van- 
ishes sufficiently quickly at infinity. It is not difficult 
to demonstrate that the same conclusion also applies to 
operator G(r,r' , E)(E — H). Since the wave functions 
— (f)ji (r) and ^^ ± (t*) are examples of functions 
which do not vanish at infinity the surface integral terms 
generated do not vanish. 

On the right hand side of Eq. (|10f) we have a purely 
scattered wave generated from cpi ? (/"). Therefore, we 
denote the result of the action of the integral operator on 
the left hand side of Eq. (jTU|l as 



dr'G(r, r'\ E ± ie)(E — H) ip±(r) 



6 (o )± , 



xi 0)± W,(13) 
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which is a part of the scattered wave xt ( r ) ■ Then for the 
part of the total scattering wave function ipk*^ developed 
from the leading term of the incident wave, 



/(0)± 



4 0)± 



(r)+xL 0)± (r), 



(14) 



we can write 

,,(0)± 



W (r) = 4> 



(r) + J dr'G(r,r';E±ie) 



(15) 



Note that although fi^^ (r) is just the leading order term 

of the incident wave, V ; i°' )± ( T ') includes scattered waves 
of all orders through the Green's function. Moreover, as 
we will see below the scattered wave in Eq. (|15|) is in fact 
the only scattered part of the total wave function ipk( r )- 
Similarly, if 4> k 15 ^ ne next-to-the-leading order term 
of the incident wave then 



(H-E^ir). 



(16) 



Therefore, for the part of the scattering wave if) k de- 



veloped from (j)^ >IC: term [i.e., for V'fe = 4>^k >± { r ) 



Xfe ± ( r )] we have 



(r) + J dr'G{r,r';E ±ie) 



xiH-E^ir') 



(17) 



By direct substitution of Eqs. ifTBf and (fT7)l into Eq. Q 

and using Eq. © we verify that ip^ ± { r ) an d V'fc 1 ^ ± ( r ) 
satisfy Eq. Obviously, continuing this procedure 

ip k (r) can ^ e formally reconstructed: 



The latter can simply be written as 



= 4>£{r)+ I dr'G(r,r , ;E±ie){H-E)<t>±(r'). 

(19) 



It is easily verified (by substitution) that this is a for- 
mal solution to Eq. QJ. Generally, this solution is not 
unique. The physical solution is the one which satisfies 
the boundary conditions specified in Eq. J5J or ©. At 
first sight the integral term in Eq. (|19fl seems to involve 
a nonintegrable singularity at the origin for the higher- 
order terms of the incident wave. However, as we will 
now demonstrate, this is not the case. 

We begin by noting as mentioned earlier that the 
leading-order asymptotic term of ip k (r) already contains 
all the scattering information we need, provided Eq. I|19|) 
is the solution which has the required asymptotic be- 
havior. In other words, for the purpose of extracting 



the scattering amplitude we have to verify that Eq. i|19|) 
asymptotically behaves as Eqs. (J2J or J3J). To this end let 
us write the full set of eigenstates of Hamiltonian H as 
{ip k , <p n }, where ip n are the eigenf unctions corresponding 
to negative discrete eigenvalues E n . Then we express the 
Green's function in its spectral decomposition as 



GirS;E±U) = fg 



3 E- k' 2 /2n±ie 
Vn{r)<p* n {r') 



E 



E-E n ± ie 

n 

Using this we can write Eq. ^19[) as 

1>tir) = 4>t{r) 

dk> {^,\H-E\4>+)^,(r) 



(20) 



(27T) 



E - k' 2 /2/i + ie 



+ 



(21) 



where dots indicate the contribution from the bound 
states. As we are interested only in the asymptotic be- 
havior when r — > oo the bound states do not contribute. 
Consider now 

(l>i\H - E\4>+) ee Jdr^,*(r)(H-E)ct>+(r), 

(22) 

the integral entering Eq. l|21|l . This is a kind of integral 
with a possible nonvanishing surface-integral component 
we mentioned earlier. Using Eq. we can write 

(?fc\H - E\4>1) = (^\-H + E + H-E\4>+) 

= -(^\Ho-Ho\K)- (23) 

Using Green's theorem this volume-integral can be trans- 
formed to a surface integral. We have 



(^\H-E\tf>+) 



fdf 



— — lim r 2 

2[l r—>oo 
1 



2^ 



lim r 2 / df 



dr 



k dr 



(24) 



Thus in the surface-integral form our integral depends 
only on asymptotic behavior of the participating func- 
tions. Now noting that 



<t>i(r) r - 
it is easy to see that 



^ 0)± (r)[l + O(l/r)], 



(25) 



fa\H - E\<f>+) = -{iP k ,\H - H \^ + ) 

= ^ k ,\H-E\4 )+ ). (26) 



4 



In combination with Eq. (|21H this result in fact means 
that 



dr'G(r, r'; E ± ie)(H - E)^(r') 

dr'G(r,r'; E ± ie)(H - E)(j> k 0)± (r'), (27) 



i.e. all the scattered wave is generated from the first term 
of the incident wave <pP . In other words, somewhat 
unexpectedly, no scattered wave is produced from the 
remaining part of (fit. This also shows that Eq. i|19|) 
does not contain nonintegrable singularities for small r' . 
Thus, using Eq. Ij26(l we can write Eq. (|21|l as 



dk' T{k',k)4> k ,{r) 
(2tt) 3 E - k ,2 /2 f i + ie 



(28) 



where dots still indicate only the contribution from the 
bound states. In the above equation we have introduced 
the notation 



T(k',k) = (^ k ,\H-E\^ )+ ), 



(29) 



anticipating that T(k', k) is the desired scattering 
T-matrix. However this requires justification that we 
now provide. In order to prove this let us first expand 
the scattering wave function ipk( r ) according to 

i>t(r) - ^V^xKfc^yfc^^f), (so) 



where 07 (fc) is the (total) phase shift. The radial func- 
tions asymptotically behave according to 

4.7T 

Xi{k,r) r ^°° — sin [kr — 7 In \2kr\ — ln/2 + cr;(fc)l . 
kr 

(31) 

Therefore, substituting Eq. I|31|) into (|3U[1 we get, after 
some algebra, 



1 — ■ x 27T 

^ (r) ~ ikr 



^ikr—i-y In \2kr 



^S(k-r) 



-ikr+i-y In \2kr 



6(k + r 



a ikr— £7 In \2kr\ 

r ik 



^( e ««^(fc) _ 1) 



I, in 



xY*(k)Y lim (?) 



(32) 



Since |2C 

27T 



_J2(e 2 * i{k) -W: m (m, m (r) = /(fc.f),(33) 



we may conclude that Eq. i1M2|l is equivalent to Eq. J5J). 
The amplitude / is still unknown, however, from here 
we get an interesting result that the first term of Eq. 



(|32|l is an asymptotic form of the distorted plane wave 
<j> k (r). Repeating the same reasoning for tpZ(r) [or us- 
ing i(j k (r) — ip^t k (r)] we establish that exactly the same 
term is the asymptotic form of 4> k °^ (r) as well. That is, 



2tt 

ikr 



^ikr—i^f In \2kr\ 



S(k 



-ikr-\-i^f In \2kr 



6{k 



(34) 



In the absence of the long-range distortion Eq. (|34|l trans- 
forms to the familiar asymptotic form of the plane wave 
(see, e.g., 0): 



2tt 

ikr 



e lkr S(k - r) 



—ikr 



5(k 



.(35) 



Thus the beauty of the result is that the notorious for- 
ward/backward singularity of the distorted plane waves 

^ , fe°' + ( T ')/ < / > i ^ _ ( r ) i s sri0wn to be no more problematic 
than a 5-function. Therefore, in the asymptotic region 
the distorted plane wave can be treated much like the 
plane wave. In this regard it is worth of mentioning that 
the Coulomb Green's function in momentum space has 
also been shown to have distorted pole singularities |2^| 
which follow from the existence of the Dollard wave op- 
erators. 

Returning to Eq. (|28() we have asymptotically, in the 
leading order, 



0i O) » + 



dk' T(k',k)^ k °)-(r) 
(2tt) 3 E-k' 2 /2n + ie 



(36) 



When r — > 00 the components involving bound states de- 
crease exponentially. Thus we have the two-body version 
of the asymptotic relationship revealed in Refs. [24l |25| . 
It states that the leading-order asymptotic term of the 
scattering wave is defined by the same (i.e., the leading- 
order) asymptotic term of the incident wave. Using 
Eq. (|34J) and evaluating the integral, taking advantage 
of the simple pole singularity of the integrand at the on- 
shell point, we have 



t/>t(r) ~ < 



/D ik-r-\-i"f In(kr-k-r) 

JL T (kr,k) — e -'7ln(2fcr). 

2tt r 



(37) 



In Eq. (|37|l we dropped the modulus sign since at this 
stage it is safe to do so. By comparing Eqs. (|3T|l and 
@ we conclude, as we set out to prove, that T(k\ k) 
introduced in Eq. (|29|l is the transition matrix (T-matrix) 
which defines the amplitude of scattering of the particles 
with initial relative momentum k in the direction of r: 



/(*• 



~2tt 



T(k',k), 



(38) 



where we used a notation k' = kr. In analogy with 
conventional scattering theory we call the new form for 
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T(k', k) as defined in Eq. i|29|) the prior form of the T- 
matrix. 

Repeating the procedure outlined above for t/>r (r) [this 
time we use the resolution of the total Green's function 
G in terms of ipf (r)] we get 

TP° st (k',k) = ($)-\H-E\iP+). (39) 

In obtaining Eq. (|39|1 we also used the reciprocity theo- 
rem, i.e. T(k',k) = T(—k,—k') 26]. Again, in analogy 
with the standard theory this new form is called the post 
form of the T-matrix. 

In deriving the result of Eq. (|26|) we saw, in particular, 
that 



T p " or (k',k) 



i^\H a -H \^ )+ ), 



1 

2^ 



lim r 2 / dr 



4 o)+ 



Or 



Or 



In a similar way we also get from Eq. Ij39|l 



T post (k',k) 



(0) 



1 

2fl 



H 
lim r 2 



,(0)-*^fe 



dr 



Or 



(40) 



.(41) 



(42) 



(43) 



Thus the scattering T-matrix conventionally given as 
volume integrals can be written equivalently in surface- 
integral forms. We emphasize that in these forms the 
T-matrix depends only on the asymptotic behavior of 
the participating functions. Therefore, knowledge of the 
scattering wave function over the entire space is not re- 
quired. In addition, the surface-integral forms are read- 
ily expanded in partial waves leading to a simple result 
containing only the limiting procedure. Therefore, these 
forms are particularly suitable for practical calculations. 
It is also interesting to note the close resemblance of these 
forms to the representation of the number of scattered 
particles crossing the surface element dr per unit time 
at large distance r. In that sense, the surface-integral 
forms further reveal the scattering amplitude as the am- 
plitude of the probability flux of particles scattered in 
direction f. 

We note that in an operator form Eq. (|19|) can be writ- 
ten as 

tffc(r) = [l + G(E±ie)(H-E)]^(r), (44) 

where G(E ± ie) = (E ±ie — H ) _1 is the Green's opera- 
tor associated with the Green's function G(r, r'; E ± ie). 
Therefore, we can introduce new generalized wave oper- 
ators according to 



fi ± = [l + G(E±ie)(H - E)]. 



(45) 



In the next section we turn to a further investigation of 
these generalized wave operators. 



III. CONSISTENCY WITH CONVENTIONAL 
RESULTS 

Our aim in this section is to show the results given 
above are consistent with conventional potential scatter- 
ing theory for short-range interactions. The existing for- 
mulation of scattering theory relies on the condition that 
interaction V(r) decreases faster than the Coulomb in- 
teraction when r — > oo [7 = in Eqs. JSJ) and ||BJ|], so 
that 

4>i(r) = ^ )± (r)^^\r) = e ik - r . (46) 

The initial unscattered wave function satisfies the 
Helmholtz equation 

(E-H Q )$\r)-< 

Then 

Mi„\ _ t/^(o). 



0. 



(47) 
(48) 



(H-E)4>r(r) = V4,%>(r), 
therefore, Eq. I|19(l takes the form 

V'fcM = < k ) (r) + J dr'G(r,r';E±ie) 
x(H-E)$\r>) 
= <f>k ] ( r ) + J dr'G(r,r';E±ie)V^ ) (r'). 

(49) 

This is the formal solution which is obtained in conven- 
tional scattering theory. In the light of Eq. Eqs. ||213|) 
and l(3*9"jl reduce to 



T*>"° r (k',k) = (?MM 0) ), 
T^(fe',fe) = (<$>\V\il>i), 



(50) 
(51) 



in agreement with the standard definitions of the T- 
matrix. Moreover, for the same reason the generalized 
wave operators introduced above reduce to the usual 
Moller (M) ones: 



= [l + G(E±ie)V]. 



(52) 

Obviously, when interaction V has a tail which does not 
disappear at infinity, the Helmholtz equation l|47|) for 

(t>k , an d consequently Eq. H48|) . are not satisfied. As a 
result Eqs. JJOJ, IHJ an d (EU are incorrect and Eq. J52J) 
is not valid for Coulomb-like potentials. 

On the other hand, when the interaction is purely 
Coulomb, we can proceed further with analytical meth- 
ods. Then we have 

(H-E)^(r)- 



[ir(kr =p k • r) 
Therefore, Eqs. (|2T?|) and transform to 



^ 0)± (r). (53) 



TP" or (k',k) = (yj k 
T post (fc',fc) = 



7 2 fc 



i° )+ >, (54) 
IV£>- (55) 



/ir(kr — k ■ r) 
7 2 fc 



fe [ir{kr + k ■ r) 
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Here ip k are known and given by 



— e ik-r -ivy/2 



r(l±i7) 
XiFi(=Fi7, 1, ±ikr — ik • r), (56) 



with \F\ being the usual confluent hypergeometric func- 
tion. At the same time Eq. I|15f) transforms to 



^ 0)± (r) = <i> 



i 0)± (r) + J dr'G(r,r';E±ie) 
7 2 fc 



jl(0)±/ ^ 



fj,r'(kr' =f fc • r') fc 



(57) 



This is to be c omp ared with the result obtained in 
Eq. (16) of Ref. [23 for the total scattering wave func- 
tion ip h : 



i 0)± (r)+ j dr'G(r,r';E±ie) 
Hr'(kr> Tk ■r>y k 1 ; ' 



(58) 



If Eq. (|58|l were true it would mean that 

dr'G(r,r';E±ie)(E-H) [^±(r') - ^ 0)± (r') 

= ^(r)-0i O)± (r),(59) 

which is, however, not correct (see Eq. 1)1 OJl and discus- 
sion following it). Based on Eq. (|58|l Barrachina and 
Macek also arrived at Eq. I|55l) . nevertheless, since the 
underlying equation was not correct this result is not 
justified in 27]. On a positive note, the matrix elements 
in Eqs. (|54|l and l|55|l (to be more precise, the complex 
conjugate of the former) have been evaluated in Ref. [27| 
in closed form. We have checked and confirm their result, 
namely 



7 2 fc 



/ir(fcr 



k ■ r 

7 2 /fc 



[ir(kr + k ■ r) 

4:Ttziz 2 r(l + ry) 
|fc'-fc| 2 r(l-i 7 ) 



4fc 2 
Ifc'-fel 2 
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(60) 



which is the well-known full on-shell Coulomb T-matrix. 
This gives additional support for the new definitions of 
the T-matrix. 

Finally, in the pure Coulomb case the generalized wave 
operators introduced earlier reduce to the wave op- 
erators obtained by Mulherin and Zinnes (MZ) pgjj : 







MZ 



l + G(E±ie) 



■y 2 k 



fj,r(kr =p k ■ r) 



(61) 



provided SV* 1 are applied to the first term of the incident 
wave ■ This clearly shows that the MZ operators 



are approximations to the corresponding full wave op- 
erators. Obviously, this approximation makes a sense 



only for asymptotically large distances where 



be- 



comes the dominant (leading-order) term. Elsewhere, the 
MZ operator cannot be relied upon. This finding may 
also explain formal problems associated with the MZ ap- 
proach [29j| . 



IV. CONCLUSION 

In this paper we have presented a generalization of 
potential scattering theory which is valid for arbitrary 
interactions including potentials with the long-range 
Coulomb tail. We obtained a new formal solution to the 
Schrodinger equation satisfying the boundary conditions 
imposed by the long-range nature of the Coulomb inter- 
action. We introduced new general definitions for the 
scattering amplitude and wave operators. We showed 
that when the interaction potential is short-ranged the 
generalized definitions of the scattering amplitude and 
wave operators transform to the conventional ones used 
in standard scattering theory. A distinctly satisfying fea- 
ture of the new forms for the scattering amplitude is that 
they do not contain divergent factors and directly give 
the physical on-shell scattering amplitude even when the 
interaction potential is long-ranged. Moreover, the gen- 
eralized wave operators are also the same for arbitrary 
potentials including the long-range interactions. There- 
fore, no modification of the theory based on renormaliza- 
tion is required. The results of the present work close the 
gap between the two different formulations of potential 
scattering theory for short- and long-range interactions. 

In conclusion we make the following comments. In the 
literature devoted to the Coulomb scattering it is custom- 
ary to separate ip^ into the 'incident' and the 'scattered' 
parts according to 



^fcW = 0*(»O + Xfc(r), 
where the incident wave is taken as [Tfj| 



(62) 



(r) = 



U^fij, 1, ±ikr — ik ■ r). (63) 



Here U is the confluent hypergeometric function of the 
second kind. The idea is driven by the fact that the con- 
fluent hypergeometric function in the regular Coulomb 
function [see Eq. i|56|) ] can be written as a sum of two ir- 
regular confluent hypergeometric functions of the second 
kind and that functions 4%(r) and Xk( r ) satisfy the first 
and the second parts of asymptotic conditions J2J or @ > 
respectively. However, (for the pure Coulomb interac- 
tion) 4>i(r) alone is a solution to the original Schrodinger 
equation, i.e.: 



(64) 



Consequently, the scattered wave xt is a solution as well. 
Thus, as a result of separation l|62|) Eq. splits into two 



7 



equations making it impossible to single out uniquely the 
important surface-integral components in the full solu- 
tion. Therefore representation (|62|l is not a satisfactory 
starting point. It also leads to other anomalies associated 
with the Coulomb problem. In particular, using Eq. I|64|) 
one can demonstrate that the Coulomb wave function is 
a solution to a homogeneous Lippmann-Schwinger equa- 
tion [30( . On the other hand the splitting (0} used in this 
work represents the logical fact that the 'unscattered' in- 
cident wave is coming from infinity and should be taken 
in a form valid at asymptotically large distances. 

We conclude by offering a few comments on the pos- 
sible uscfulcncss of our methods for three-body systems 
with long-range interactions. Rigorous scattering theory 
for a system of three particles valid for short-range poten- 
tials was given by Faddeev |3ll |32| ■ For the charged par- 
ticles with the long-range Coulomb interaction the the- 
ory has faced apparently insurmountable difficulties. The 
problem is that the Faddeev equations are not compact 
in the presence of Coulomb interactions. In other words 
these equations cannot be solved using standard numeri- 
cal procedures. Some progress has been made in dealing 
with aspects of the three-body problem with Coulomb- 
like potentials. In particular, a renormalization method 
based on screening 0, Q has been implemented success- 
fully for the case when two particles are charged [3^, OSI • 
The method has also been extended to two-fragment re- 
actions in a system of three charged particles HHl • 
Though Dollard's time-dependent approach 0, @] is for- 
mally valid for arbitrary multichannel collisions includ- 
ing three-body problem, no practical time-independent 



renormalization method exists that is valid for a system 
of three charged particles above the breakup threshold. 
The problem is that above the threshold the Coulomb 
three-body system possesses essentially different types 
of singularities and the two-particle renormalization pro- 
cedures are not sufficient to guarantee compactness of 
the equations 0, |3(| |37|]. Thus, on one hand there are 
no integral equations yet known for collisions of more 
than two charged particles that are satisfactory above the 
breakup threshold an d on the other hand, there is 
also neither theoretical proof nor practical evidence that 
renormalization approach can be applied to the Faddeev 
equations for the genuine three-body Coulomb problem. 
This is a rather disturbing situation especially for the 
atomic three-body problem where all three particles are 
charged. Therefore, generally speaking, it would be use- 
ful to formulate three-body scattering theory in a manner 
that does not require renormalization so that the afore- 
mentioned modifications are, in a certain sense, unnec- 
essary. We are confident that the method proposed here 
for solving the two-body problem that was free of the 
usual Coulomb anomalies can be profitably applied to 
the proper formulation of the three-body rearrangement 
theory for long-range interactions. 
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